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E865 at the Brookhaven National Laboratory AGS collected about 70,000 K^ 3 events to measure 



the K2 



branching ratio relative to the observed K + 



7T", K + — > 71° /J, 1 V 

decays. The 7r u in all the decays was detected using the e + e~ pair from 7r° — * e 
photons were required. Using the Particle Data Group branching ratios Q for the normalization 
decays we obtain BR(K^.<) = (5.13 ± 0.02 sta t ± 0.09 a!/ s ± 0.04 norm )%, where K^, » includes the 
effect of virtual and real photons. This result is w 2.3a higher than the current Particle Data Group 
value. Implications for the V us element of the CKM matrix, and the matrix's unitarity are discussed. 

PACS numbers: 13.20.Eb, 12.15.Hh 
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The experimentally determined Cabibbo-Kobayashi- 
Maskawa (CKM) matrix describes quark mixing in the 
Standard Model framework. Any deviation from the ma- 
trix's unitarity would undermine the validity of the Stan- 
dard Model. One unitarity condition involves the first 
row elements: 
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where a non-zero value of S indicates a deviation from 
unitarity. The V u d element is obtained from nuclear and 
neutron decays. Kb, from the semileptonic decays of 
B mesons is too small to affect Eq. ^ The V us 
element can be determined either from hyperon, K — » 
■k^lv(K^) or from K — > irev (K e s) decays. However,/^ 
decays provide a smaller theoretical uncertainty 0, 0, 
The most precise value of V u d obtained from the nuclear 
superallowed Fermi beta decays leads to S = (3.2 ± 1.4) • 
10 -3 0, a 2.3cr deviation from unitarity. 

Both experimental and theoretical efforts to improve 
the determination of Kd continue. Theoretical contri- 
butions to V us were reevaluated recently 0, 0, 0, Q, but 
there has been little new experimental input on the K ij 
branching ratio. Since the V 2 d and V 2 S uncertainties are 
comparable, a high statistics measurement of the K* 3 
branching ratio (B.R.) with good control of systematic 
errors is useful. 

The bare (without QED corrections) K^ 3 decay rate 
mill! is: 

dT(K+) = C7(t)|K s | 2 |/ + (0)| 2 [l + K^Ydt (2) 

where t = {Pk — P-rr) 2 , C(t) is a known kinematic 
function, and /+(0) is the vector form factor value at 



t = 0, determined theoretically 0, 01 Two recent 
experiments [ltl IbH ] give A+ (the form factor slope) mea- 
surements consistent with each other and with previous 
measurements. An omitted negligible term contributing 
to Eq. [21 contains the form factor /_, and is proportional 
to Mil Ml. 
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FIG. 1: Plan view of the E865 detector with a simulated 
K + — > Tv°e + u decay followed by it — » e + e~-f. 



E865 .12| searched for the lepton flavor violating de- 
cay K + — ► 7r + /i + e~. The detector (Figure resided in 
a 6 GeV/c positive beam [l2^. For the K+ running, the 
intensity was reduced by a factor of 10, to 10 7 kaons, 
2 x 10 8 protons, and 2 x 10 8 7r per 2.8 second spill. The 
beam was intentionally debunched at extraction to re- 
move rf structure at the experiment. The first dipole 
magnet separated particles by charge, while the second 
magnet together with four multiwire proportional cham- 
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bers (MWPCs: P1-P4) formed the spectrometer. The 
particle identification used the threshold multichannel 
Cerenkov detectors (CI and C2, each separated into left 
and right volumes, for four independent counters) filled 
with gaseous methane (Cerenkov threshold jt ~ 30 and 
electron detection efficiency e e w 0.98 0), an electro- 
magnetic calorimeter an( i a rnuon detector (not used 
for the K^ 3 measurement). The D and A scintillator ho- 
doscopes gave left/right and crude vertical position. 

The 7T° from the kaon decays was detected through 
the e + e~ from the ir° — > e + e _ 7 decay, with the 7 
detected in some cases. To eliminate the uncertainty 
(2.7%) of the 7T° — > e+e _ 7 B.R., and to reduce system- 
atic uncertainty we used the other three major decay 
modes with a tt° in the final state (K + — > 7r + 7r°(if 7 J 2 ), 
K^ 3 , K + — » 7r + 7r 7r (A"^3)) for the normalization sample 
("Kdal"). 

The if^ data was collected in a one-week dedicated 
run in 1998, with special on-line trigger logic. 

The Kdal and K^ 3 data were collected by the "e + e~" 
trigger, which was designed to detect e + e~ pairs and re- 
quired at least one D-counter scintillator slat on each 
(left and right) side of the detector and signals from each 
of the four Cerenkov counters. The Cerenkov efficiency 
trigger required only 3 out of 4 Cerenkov counters (no D- 
counter requirement). The "TAU" trigger, requiring only 
two D-counter scintillator hits (one left, and one right), 
collected events for the K + —> it + it + it~ (K t ) sample, to 
study the detector unbiased by Cerenkov requirements. 
About 50 million triggers were accumulated, ~ 37 mil- 
lion in the "e + e~" trigger. About 75% of "e + e _ " triggers 
included accidental tracks, often a \i from high momen- 
tum K — * [iv or 7r — ► \w decays partially satisfying the 
Cerenkov requirement. 

Off-line reconstruction used the spectrometer only. 
The Cerenkov and D counter efficiencies were obtained 
from the Cerenkov efficiency triggers. The redundancy 
of the MWPCs (4 planes/chamber) and track reconstruc- 
tion was used to extract MWPC efficiencies. The absence 
of the electromagnetic calorimeter from the trigger al- 
lowed its efficiency determination. Each efficiency was 
measured over its relevant phas e sp ace. 

Relevant kaon decay chains [13| were simulated with 
GEANT [lj] (including decays of secondary pions and 
muons). For A+ = 0.0278 ± 0.00190 was used. 

The radiative corrections to the K^ 3 decay phase-space 
density || were used. The K+ 3 (inner bremsstrahlung) 
decays outside the K^ 3 Dalitz plot boundary were explic- 
itly simulated [tj. For ir° —> e + e _ 7 decay, radiative cor- 
rections were taken into account according to Ref. [l5| . 
Measured efficiencies were applied 13], and accidental de- 
tector hits (from reconstructed K T events) were added. 
About 10% of both the K^ 3 and Kdal samples had extra 



decay volume (no requirement for exactly three recon- 
structed tracks was applied) , for the three tracks to cross 
the active parts of the detector, for the low (M ee < 0.05 
GeV) mass e + e~ pair to be identified in the Cerenkov 
counters, and for the second positive track to have less 
than 3.4 GeV/c momentum. The momentum cut re- 
jects events where /i + or tt + from Kdal decays is above 
Cerenkov threshold and can be identified as e + . A ge- 
ometric Cerenkov ambiguity cut rejected events (27%, 
15%, 25%, and 35% for K+ 3 , K+ 2 , K+ 3 , and K+ 3 respec- 
tively) where the Cerenkov counter response could not 
be unambiguously assigned to separate tracks[l3^. 

The K^ 3 sample was then selected by requiring the 
second positive track to be identified as e + in 2 of the 
3 electron detectors: CI, C2, or the calorimeter, each 
with e e w 98%. Events entering the Kdal sample had no 
response in at least one of the two Cerenkov counters. 
These criteria minimized systematic uncertainties |l3j . 
but resulted in a small overlap, w 3% of the K^ 3 sample 
and « 0.3% of the Kdal which was accounted for in the 
B.R. calculation. The Ki, acceptance is « 1.2%. The 
K^ 3 acceptance « 0.7% |13|. somewhat lower because 
of the lower average e + momentum in the K^ 3 decay. 
The overall acceptance level of 1% can be approximately 
understood by assuming a factor of three loss for each 
charged particle, 30 % for the Cerenkov ambiguity, and 
approximately a factor of 2 for other cuts. Final accep- 
tances for the three modes in the Kdal sample differed by 
< 4% taking into account that either of the 7r°s from K^ 3 
can decay into e + e~"f. The final K^ 3 and Kdal samples 
were 71,204 and 558,186, respectively. Figure El shows 
some relevant spatial distributions. 




reconstructed tracks. 

Selection criteria, common to K 



and Kdal, included 
requirements for a good quality three track vertex in the 



< 3 



FIG. 2: Distributions of X and Y positions of the first positive 
track (not e + from the 7r° decay) for the selected K^ 3 and 
Kdal samples. X and Y positions are measured at the end of 
the second pair of the Cerenkov counters (C2). Histograms 
represent Monte Carlo; points with errors represent data. 

Contamination of the K^ 3 sample by other K + de- 
cays occurred when n + or fi + from Kdal decays were 
misidentified as e + , or as a result of n° — > e + e~e + e~ . 
Contamination due to secondary particle decays was es- 
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timated to be at the level of 0.1%. About 8% of final state 
pions decayed into muons inside the spectrometer. The 
careful MWPC simulation gave good agreement of recon- 
structed track x 2 and vertex distributions between data 
and Monte Carlo. No tight track x 2 cuts were applied, 
and the systematic uncertainties estimated by variation 
of the vertex cuts were included in the final result. The 
check of B.R.( AT T /Kdal), described below, also tests the 
final state it and p, decays. 

Total contamination of the Ke3 sample was estimated 
to be (2.49 ± 0.05 s t at ± 0.32 sys )%, with the systematic 
uncertainty caused by the simulation accuracy of the CI 
and C2 response to 7r + and fi + . Contamination due to 
overlapping events was (0.25±0.07)% and (0.12±0.05)% 
of the Kdal and K+ 3 respectively. Figure shows the en- 
ergy distribution in the calorimeter from the e + in the 
-ftfij sample. The contamination is manifest in the min- 
imum ionization spike at 250 MeV. The small excess of 
data in the spike agrees with our contamination uncer- 
tainty estimate. 

The final K^ 3 sample included «30% of events with a 
fully reconstructed 7r°s. We used the 7r° information as a 
consistency check. Not requiring 7r°s in our main analysis 
minimized the uncertainty arising from photon detection 
and reconstruction in the calorimeter, but increased vul- 
nerability to contamination from upstream decays and 
photon conversion. Upstream decays whose photon pro- 
duced pairs before the decay volume (evacuated to about 
10~ 8 nuclear interaction length) were suppressed by re- 
quiring the three track vertex to be more than two meters 
downstream of the decay volume entrance. In addition, 
the results obtained from the two independent samples, 
one with and one without the ir° reconstructed, did not 
show a statistically significant discrepancy. 

The K^ 3 statistical precision is 0.4%. The systematic 
error estimate, summarized in Table [U was determined 
from the B.R. stability under variation of reconstruction 
procedure, selection criteria, assumed detector efficien- 
cies, and subdivision of both iCt. and Kdal samples |l3|. 
No significant correlations between any of the different 
systematic uncertainties were observed. 

The two largest contributions to the systematic error 
come from the discrepancies 0] between data and Monte 
Carlo in the momentum (Figure and spatial distribu- 
tions. These errors were determined by dividing the K^ 3 
and Kdal events into two roughly equal subsamples, using 
the relevant parameters, and observing the variation of 
the result [13J. The errors were found to be uncorrelated. 
The sensitivity of the vertical spatial discrepancy to the 
MWPC alignment and of the momentum discrepancy to 
the spectrometer parameters is indicative of their pos- 
sible origins 0| . The Z-vertex position is also sensitive 
to the magnetic field, but has a smaller systematic er- 
ror contribution as determined from both upstream and 
downstream cuts in Z. 

As an additional consistency check, we estimated the 
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FIG. 3: Energy deposited in the calorimeter by the second 
positive track from the selected K^ 3 sample (e + which is not 
from the low mass e + e~ pair). No calorimeter information 
was used for the e + identification. Markers with errors repre- 
sent data; the histogram is simulation. 



TABLE I: Systematic uncertainty sources and estimates of 
their respective contributions to the final result's uncertainty. 
The total systematic error is the sum (in quadrature) of the 
individual contributions. 



Source of systematic error 


Error 
estimate 


Magnetic field uncertainty 


0.3% 


Vertex finding and quality cut 


0.6% 


Vertex position cut 


0.4% 


Cerenkov Ambiguity Cut 


0.3% 


M ee Cut 


0.2% 


Detector Aperture 


0.2% 


(Tv/fi) + identification 


0.04% 


MWPC efficiencies 


0.2% 


D counter efficiencies 


0.15% 


Cerenkov efficiencies 


0.3% 


Contamination of the selected samples 


0.3% 


Removal of extra tracks 


0.2% 


Vertical spatial/angle distributions discrepancy 


0.8% 


e + /e~ momentum distributions discrepancy 


1.3% 


K^ 3 trigger efficiency 


0.1% 


Uncertainty in the K^ z form factor slope 


0.1% 


Total error 


1.8% 



AV/Kclal B.R.. The result was (1.01 ± 0.02) xthe PDG 
ratioQ, (the theoretical prediction 0] was used for the 
7T° — » e + e~7 decay rate). The 2% error was dominated 
by the uncertainty in the prescale factor of the TAU trig- 
ger. A second consistency check compared the K^ 3 B.R. 
from 1998 and 1997 data. The 1997 K+ data used a trig- 
ger that required calorimeter hits, and A and D-counters. 
That trigger neither allowed measurement of these detec- 
tor efficiencies nor of the trigger efficiency. While we did 
not use the 1997 data for our final result, the 1997 K^ 3 
branching ratio was statistically consistent (within one 
sigma) with that from 1998. This agreement is important 
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FIG. 4: Reconstructed momentum of the e from the low 
mass e + e~ pair from the selected K^ 3 and Kdal samples. His- 
tograms represent Monte Carlo; points with errors represent 
data. Plots on the right show the bin by bin Monte Carlo to 
data ratio. 



since the momentum spectrum discrepancy between data 
and Monte Carlo in the 1997 data is qualitatively differ- 
ent from 1998|L3j. A preliminary reconstruction version 
was used for the 1997 data, without the final magnetic 
field and detector alignment. This bolsters our intuition 
that the discrepancies in decay product momenta and 
spatial distributions, which dominate the systematic un- 
certainties, reflect our imperfect knowledge of the mag- 
netic field and detector positions but do not bias our 
result beyond our estimated systematic errors. 

We estimated the form factor slope A+ from both 1998 
and 1997 K+ data l.T. We obtained: A+ = 0.0324 ± 
0.0044 stQt for 1998, and A+ = 0.0290 ± 0.0044 stat for the 
1997 data, both consistent with the current PDG fit. 

After contamination subtraction |l.3j| . our result 



is 



BR(K+ h) )/(BR(K+ 2 ) + BR(K+ 3 ) + BR(K+ 3 )) = 
0.1962 ± 0.0008 s4at ± 0.0035 sys , where K+ h) includes 
all QED contributions (loops and inner bremsstrahlung). 

3 + e" pair 

from 7T° 



As noted above, the n° was detected using the e^e' 
'2, and n o photons were required. 
Using currentyj Kdal B.R.'s we infer BR(K + 



(5.13±0.02 fltat ±0.09 sya ±0.04 r , 



where the normal- 



ization error was determined by the PDG estimate of the 
Kdal B.R. uncertainties. This result does not include the 
correction due to the correlation of the PDG kaon decay 
ratios, since it was estimated to be small compared to the 
systematic error. The PDG fit to the previous K + decay 
experiments yields BR(K+ n®e+v) = (4.87 ± 0.06)% 
1], ~ 2.3c lower than our result. 

Radiative corrections for decays inside the K^ 3 Dalitz 
plot boundary were estimated to be —1.3% using the 
procedure of Ref. 0; K^-y decays outside the Dalitz 
plot boundary gave +0.5%. Thus the total radiative 
correction was —0.8% resulting in the bare BR(K^ 



^3) 



(5.17±0.02 stat ±0.09 sys ±0.04 r , 



hancement factor S E w(M p ,M z ) = 1.02320113, and 
our result for the bare K^ 3 rate we obtain \V us f+(0)\ = 
0.2243 ± 0.0022 rate ± 0.0007 A+ , which gives = 
0.2272 ± 0.0023^^ 0.0007 A+ ± 0.0018 /+(0) if /+(0) = 
0.9874 ± 0.0084011. With this value of V us and V ud 
from superallowed nuclear Fermi beta decays^, 8 — 
0.0003 ±0.0016. 



This result is consistent with CKM unitarity, but 
increases the discrepancy with the V us from K® 3 de- 
cay if extracted under conventional theoretical assump- 
tions about symmetry breaking. K e ^ measurements in 
progress (CMD2, NA48, KLOE)0 should help to clarify 
the experimental situation. 
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Using the PDG value for Gf, the short-distance en- 
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